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5 ■ Abstract 

I , We discuss a many-body Hamiltonian with two- and three-body interactions 

^ , in two dimensions introduced recently by Murthy, Bhaduri and Sen. Apart 

^ from an analysis of some exact solutions in the many-body system, we ana- 

^ , lyze in detail the two-body problem which is completely solvable. We show 

H : 

that the solution of the two-body problem reduces to solving a known differ- 
ential equation due to Heun. We show that the two-body spectrum becomes 
remarkably simple for large interaction strength and the level structure resem- 
bles that of the Landau Levels. We also clarify the "ultraviolet" regularization 
which is needed to define an inverse-square potential properly and discuss its 
implications for our model. 
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I. INTRODUCTION 



Recently, there has been a revival of interest in the area of exactly solvable models in one 
and higher dimensions. A celebrated example of a solvable many-body system is the well- 
known Calogero-Sutherland model (CSM) in one dimension The model has found wide 
application in areas as diverse as quantum chaos and Fractional Statistics. The particles 
in the CSM are confined in a one-body oscillator potential or on the rim of a circle, and 
interact with each other through a two-body potential which varies as the inverse-square of 
the distance between particles. The CSM and its variants in one dimension, like the Haldane- 
Shastry model for spin chains , have provided the paradigms to analyze more complicated 
interacting systems. A characteristic feature of the CSM is the structure of the highly 
correlated wave function. The correlations are built into the exact wave function through a 
Jastrow factor ( I" for any pair of particles denoted by The exponents 
on the correlator are related to the strength of the inverse-square interaction. Notice that 
this factor is antisymmetric (symmetric) in particle labels for A = 1(0) and vanishes as the 
two particles approach each other. A generalization of this in two dimensions is to be found 
in Laughlin's trial wave function |P where the correlations are built in through the factor 
{zi — Zj), where Zi are the particle coordinates in complex notation. The corresponding 
Hamiltonian for which the Laughlin wave function is an exact ground state has not been 
analyzed to the same degree of detail as the CSM. It is known that it is the ground state for 
a Hamiltonian describing spin polarized electrons in the lowest Landau level with a short- 
range repulsive interaction P|. It is also known that such correlations are present in the 
exact ground state of a spin Hamiltonian in two- dimensions. The anyon Hamiltonian p| 
in two-space dimensions is another example where the Jastrow correlation appears . While 
the two-anyon problem is exactly solvable, the many-anyon problem is not. For a system 
of anyons confined in an oscillator potential many exact solutions and their properties are 
known but, unlike the CSM in one- dimension, the analytical solution of the full many-body 
problem is not tractable [0. It is therefore of great interest to find models analogous to 
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the CSM in higher dimensions. 

In a recent paper ||Tl|, three of us proposed a model in two-space dimensions with non- 
trivial two- and three-body interactions which could be solved exactly for the ground states 
and some excited states. It betrayed some similarity to both CSM in one-dimension and 
the anyonic model in two-dimensions through the spectrum. The model was devised by 
noting that in two dimensions there exists another form of the pair correlator with which a 
Jastrow-type many-body wave function may be constructed, namely 

Xij = XiHj - XjHi . (1) 

The correlation is by definition antisymmetric and goes to zero as two particles approach 
each other. In addition, it introduces zeros in the wave function whenever the relative angle 
between the two particles goes to zero or vr. The difference with the Jastrow-Laughlin form 
is also significant; Xij in (|l]) is a pseudo-scalar. Unlike the Laughlin type of correlation, 
it does not impart any angular momentum to the two dimensional wave function . One 
important drawback of this correlation is that it is not translationally invariant unless the 
radial degrees of freedom is frozen. The model Hamiltonian has solutions which have this 
correlation built in. Intuitively the correlation can be understood easily by imagining objects 
with associated "arrows". The arrows cannot be oriented either parallel or anti-parallel to 
each other. The model has some interesting features and it would be of great interest to 
find physical systems which incorporate these features. 

In this paper we elaborate on our earlier results |jll| and present several new results. In 



Sect. II, we discuss the many-body Hamiltonian and display some of the exact solutions 
and their structure. The similarities between the spectrum of these exact solutions and the 
spectrum of CSM is quite remarkable. Further, when projected on to a circle the model 
reduces to a variant of the trigonometric Sutherland model. In this limit the model also 
has translational invariance. In Sect. Ill, we discuss the two-body problem in detail and 
show that the solutions of the two-body problem are described by the Heun equation. In 
particular, the spectrum becomes very simple for large values of the interaction strength. 



The singular interaction discussed in this paper requires a careful treatment in the region 
near Xij = 0; this is discussed in the Appendix |jl2[. Sect. IV contains a discussion and 
summary. 



II. MANY-BODY HAMILTONIAN AND SOME EXACT SOLUTIONS 

For the sake of completeness we recall first the Hamiltonian and some of its properties 
proposed earlier [jlT|. We also clarify some points which were not made explicitly clear in 
the earlier paper. The iV-particle Hamiltonian which we consider is given by 

1=1 1=1 i.j ■>',J,k 'J ''^ 

where X^j is given by ([1|); gi and g2 are dimensionless coupling strengths of the two- and 
three-body interactions respectively. While gi and g2 can be independent of each other in 
general, for the type of solutions involving the correlator in (||) they are not. We will specify 
their relationship shortly. The particles are confined in a one-body oscillator confinement 
potential. The Hamiltonian is rotationally invariant and manifestly symmetric in all particle 
indices. As in the CSM, we may scale away the mass m and oscillator frequency u by scaling 



all distances — >■ ^Jmui/h fi, and measuring the energy in units of huo. This is done by 
setting h = m = u = 1. In these units the Hamiltonian is given by 

T N -1 N n ^ n N ^ ^ 

= -5Ev| + iEff + f Ej^ + f E f^- (3) 

^ i = l ^ i=l ^ ^.3 ^ij ^ ^,J,k ^ij^ik 

Note that the total angular momentum operator 

L = Yl i^iPy^ - ViP^.) (4) 

i 

commutes with the Hamiltonian since it is rotationally invariant, and may therefore be used 
to label the states. The Hamiltonian is invariant under parity x —x and y y. In 
addition, for any i, the Hamiltonian is invariant under the transformation —fi and 

"f^k — >■ Tk for all k ^ i. This D2N invariance is special to this system, and we are not aware of 



any other interacting many-body Hamiltonian which has this symmetry. The consequences 
of this will be discussed explicitly in the two-body problem where this is related to the 
supersymmetric properties of the system. 

We will consider both bosonic and fermionic systems governed by the Hamiltonian 
i.e., wave functions which are totally symmetric and antisymmetric respectively. It will turn 
out that certain calculations (for example, in the two-body problem) simplify if we do not 
impose any symmetry to begin with. 



A. The exact bosonic ground state 

We first obtain the exact bosonic ground state of this Hamiltonian. As an ansatz for the 
ground state wave function, consider a solution of the form 

N 1 Af 

Mx^,y^) = n exp (--E^l') • (5) 

i<j 1=1 

Clearly \E'o correctly incorporates the behavior of the wave function in the asymptotic re- 
gion \fi\ — s> oo, and \E'o is regular for g > 0. In general we insist that our solutions have 
this asymptotic form; the conditions under which this is valid will be specified later. The 
eigenvalue equation now takes the form 

1 ^ I ^ f -fu 

H^o = Li9i - 9(9 - I)) + 0^92 - g') E -^^ + 9N{N-l) + N]^o. (6) 

Therefore \l/o is the exact many-body ground state for an arbitrary number of particles of 
the Hamiltonian if 

9i = 9(9 - 1) and g2 = g^ . (7) 

Since (7 > 0, we have gx > —1/4 and g2 is positive definite. Note that the range of gi is 
identical to the one obtained in the CSM. The ground state energy is now given by 

Eo = N + gN{N - 1) . (8) 

Note that this has exactly the form of the ground state energy of the CSM. 



Since g determines both gi and g2 uniquely, we will regard g as the fundamental param- 
eter of the Hamiltonian which determines the strength of the interaction. In other words, 
we demand that the ground state should be given by (^), and we define the Hamiltonian to 
ensure that. It turns out that such a definition requires some special care in the vicinity of 
Xij = (called the ultraviolet region below). The Appendix will discuss this for the case of 
two particles. Note that in the two-particle case, the Hamiltonian only contains the param- 
eter gi = g{g — 1) and not g2- As a result, for every value of gi in the range —1/4 < gi < 0, 
the bosonic ground state energy as given hj Eq = 2 + 2g has two possible values; these two 
possibilities correspond to different potentials in the ultraviolet region. This is somewhat 
unusual but it is not uncommon for singular potentials. The same thing also happens in 



the CSM even for the A^-body problem; see for example Ref. [|T3[. We discuss this issue in 
detail in the Appendix where we show that the ultraviolet regularization is determined by 
the parameter g rather than by gi. 

We emphasize that our objective here is not to find the general solutions for arbitrary 
gi and g2, but to find a Hamiltonian whose solutions have the novel correlation in Eq. (|Tp 
built in. In general, if gi and g2 are independent, the Hamiltonian will have a ground state 
different from the one given above. Our procedure is therefore similar to the many-anyon 
problem where also there are two- and three-body interactions, but the strengths are related 
to a single parameter. With the form of gi and g2 given in (^, the solution found above is 
indeed the lowest energy state. 

A neat way of proving that we have indeed obtained the ground state can be given using 



the method of operators [14|. To this end, define the operators 



Qxi = - ixi + i g 51 ^ ' 
QU = Px, + ixi - i g I] ^ , 

Ql = Py. + ^y^ + ^9 E 1^ • (9) 



It is easy to see that the Q's annihilate the ground state in Eq. 

Q^^^o = and Qy^^o = 0. 
The Hamiltonian can now be recast in terms of these operators as 

I E [QlQ^.+QlQy] = H ^ Eo, (10) 

i 

where Eq is given by Eq. (H). Clearly the operator on the left hand side is positive definite 
and annihilates the ground state wave function given by Eq. @. Therefore Eq must be the 
minimum energy that an eigenstate can have. 

As we remarked earlier, the ground state of the Hamiltonian is bosonic. The ground 
state of the Hamiltonian for a fermionic system is not easy to determine analytically (for 
g > 0). The problem here is analogous to a similar problem in the many-anyon Hamiltonian 
|15| , p!6[| . In Sect. HI, we will determine the fermionic ground state energy for two particles 
both numerically and to first order in g using perturbation theory near g = and show that 
it has quite unusual behavior. 



B. Spectrum of excited states 

While we have not been able to find the complete excited state spectrum of the model, 
the eigenvalue equation for a general excited state may be obtained as follows. From the 
asymptotic properties of the solutions of the Hamiltonian in Eq. (^, it is clear that \E' has 
the general structure 

*(xi,2/i) = '^o{x^,yi) <!>{xi,yi) , (11) 

where \E'o is the ground state wave function. Obviously if $ is a constant we recover the 
ground state. In general $ satisfies the eigenvalue equation 

[ - i E V? + Ef. . V. + , E - A) i,^iE-E.),. (12) 
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It is interesting to note that while gi is zero both at g = and 1, the term containing g in 
the above expression is zero only when g = 0. This is so because of the boundary condition 
that the wave functions must vanish as \Xij\^ for nonzero g. 

We first discuss the exact solutions of the above differential equation. This is easily done 
by defining the complex coordinates 

z = X + iy and z* = x — iy , (13) 

and their partial derivatives 

d = d/dz = ^ {d/dx—id/dy) , 
d* = d/dz* = ^ {d/dx + id/dy) . (14) 

In these coordinates, the differential equation for $ reduces to i7$ = {E — Eq)^, where 

H = -2^28,8* + EM. + ^^D + 2^? E ^4^^ . (15) 

In addition, $ is an eigenstate of the total angular momentum operator, L$ = /$. We can 
now classify some exact solution according to their angular momentum, 
(a) / = solutions: Define an auxiliary parameter 



t = E 



ZiZ^ , 



and let $ = $(t). This has zero total angular momentum. The differential equation for $ 
reduces to 

where b = Eq and a = {Eq — E)/2; Eq is the energy of the ground state. The allowed 



solutions are the regular confluent hypergeometric functions [17 



<l>(t) = M(a,6,t) . (17) 
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Normalizability imposes the restriction a = —rir, where is a positive integer; then 

is a polynomial of degree (the subscript 'r' denotes radial excitations as discussed later). 

The corresponding eigenvalues are 

E = Eo + 2nr . (18) 



This class of solutions was discussed before in [p |. 

(b) Z > solutions: Let 

i 

and let $ = $(tz). The total angular momentum is not zero. All the mixed derivative terms 
in Eq. (|15]) drop out, and we get the differential equation 

2t.^ = (E-Eo)^. (19) 

This is the well known Euler equation whose solutions are just monomials in t^. The solution 
is given by 

HQ = C , (20) 

and the total angular momentum is / = 2m. The eigenvalues are 

E = Eo + 2m = Eo + l . (21) 

(c) / < exact solutions: Let 

i 

and let $ = ^(t^*). Once again the differential equation for $ reduces to 

= {E~Eo)^. (22) 
This is similar to the previous case. The solution is given by 

$(t.O = , (23) 

and the total angular momentum is / = —2m. The eigenvalues are 

10 



E = Eo + 2m = Eo-l . (24) 

(d) Tower of excited states: One can now combine solutions of a given / in cases (b) or (c) 
with the solutions in (a), and get a new class of excited states. Let us define 

= $i(t)$2(t.) , (25) 

where $1 is the solution with / = 0, $2 is the solution with / > 0, and t and tz have been 
defined before. The differential equation for $ is again a confluent hypergeometric equation 
given by 

t— + (i-t)--ai = 0, (26) 

where b = Eq + 2m and a = {Eq + 2m ~ E)/2. The energy eigenvalues are then given by 

E = Eo + 2nr + 2m = Eq + 2nr + I . (27) 

One may repeat the procedure to obtain exact solutions for a tower of excited states with 
/ < solutions. As we shall see below, the existence of the tower is a general result applicable 
to all excited states of which the exact solutions shown above form a subset. We notice that 
these solutions bear a remarkable resemblance to the many-anyon system where a similar 



structure exists for the known class of exact solutions p!0|| . 

(e) A general class of excited states: One can combine the solutions of all the three classes 
(a), (b) and (c) to obtain an even more general class of solutions. Consider the polynomial 

P(ni,n2,n3) = C*' , (28) 

where the rii are non-negative integers. Using the form in ([T5|), one can show that 



HP{ni,n2-,rfi) = 2 (ni + n2 + n3)P(ni, n2, ns) - 8n2n3P(?2i + 1, n2 - 1, ns - 1) 

- 2ni[ni + 2n2 + 2n3 + gN{N -l)]P{ni-l,n2,ns) . (29) 

Using this one can show that there is an exact polynomial solution, whose highest degree 
term is P{ni, n2, n^). The energy of this solution is 
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E = Eo + 2(ni + n2 + n3) , (30) 

and the angular momentum is / = 2{n2 — ns). 

While there may be more exact solutions, we do not know of a simple way of solving for 
them. We can however glean some general features as follows. The coordinates {xi,yi) can 
be separated into one 'radial' coordinate t = as above and 2N — 1 'angular' coordinates 

collectively denoted by fij(say). Then, the Eq. ( [T^ ) can be expressed as 

9^$ 9$ 1 1 

+ (Eo-t)— - ^^<f + 2(^-^0)$ = 0, (31) 

where C = T>2+g'Di, and P„ is an n^'^-order differential operator which only acts on functions 
of the angles Vti. In particular, X>2 is the Laplacian on a sphere of dimension 2N — 1. Next 
we note that the $ can be factorized in the form 

$(a;,,|/,) = R{t) r(a) , (32) 

where F, generalized spherical harmonic defined on the 2N — 1 dimensional sphere S*^^"^, 
satisfies the eigenvalue equation CY = \Y. (This is the hard part of the spectral problem, 
to find the eigenvalues A). We now define 



/i = ^{Eo-iy + AgX - (Eo-l). (33) 

Further if we write R{t) = t^^'^R(t), then R satisfies a confluent hypergeometric equation 

(PR /, ^ dR „ /„,\ 

i_ + (t_i)__,fl = 0. (34) 

where b = Eo + fi and a = {Eq + fi — E)/2. The admissible solutions are the regular confluent 
hypergeometric functions, R{t) = M{a,b,t). Normalizability imposes the restriction a = 
—rir, where is a positive integer. Then R{t) is a polynomial of degree n^, and it has 
nodes. The energy of this state is given hj E = Eo + ^ + 2nr. We see that for a given value of 
/i, there is an infinite tower of energy eigenvalues separated by a spacing of 2. As remarked 
earlier, this is reminiscent of what happens in the case of anyons. The tower structure and 
the angular momentum are useful in organizing a numerical or analytical study of the energy 
spectrum. Since the radial quantum number and the angular momentum / are integers, 
they cannot change as the parameter g is varied continuously. 
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C. Relation to Sutherland model 



It may be of interest to note that the model reduces to a variant of the Sutherland model 
in one dimension. In this limit, therefore the model is exactly solvable. Restricting the 
particles to move along the perimeter of a unit circle in the Hamiltonian (|^) without the 
confinement potential, we get 

-jEa^ + f + | E |l+cot(»,-»,)cot(»,-ft)]. (36) 

1=1 ^ X.J V t J/ %.j.k 

since Xij = — sm{9i — 9j) now. Using the identity 

f: cot(0.-^^,)cot(^^.-^^.) = - ^ (36) 

i.j.k 

we immediately recover an analog of the trigonometric Sutherland model, but shifted by the 
constant g2 N{N — 1)(A^ — 2)/3. Note, however, that the potential in ( pS] ) depends on the 
function sin(^^j — dj), rather than the chord-length which is proportional to sin[(^^j — 0j)/2] . 
Interestingly the wave function has twice the periodicity of the Sutherland model solutions- 
the wave function vanishes whenever the particles are at diametrically opposite points on a 
circle or at the same point. 

III. THE TWO-BODY PROBLEM: COMPLETE SOLUTION 

While we have not been able to solve the many-body problem completely, the two- 
body problem in our model is exactly solvable. We demonstrate this by going over to the 
hyperspherical formalism first proposed in two dimensions by Kilpatrick and Larsen [|18| (see 
also |]15[). We discuss some of the properties of the two-body spectrum. We also explicitly 
show that the two-body problem is integrable. It is important to note that the two-particle 
interaction is sufficiently singular that a careful treatment is required in order to define the 
problem completely consistently; this is described in the Appendix. 

The two-body Hamiltonian is given by 
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H = -^[V? + V^] + I[r1 + r1] + |^ , (37) 

where X = Xiy2 — X2yi. The two-body problem is best solved in the hyperspherical coor- 
dinate system which allows a parameterization of the coordinates ri,r2 in terms of three 
angles and one length, (i?, 9, 0, ■0) as follows: 

^1 +Wi — R (cos 6 cos (f) — isva.6 sin 0) exp(i'0) , 

X2 + W2 = R {cos9 sm{(f)) + i sin 9 cos cf)) cxp{itp) . (38) 

We may regard {R, 9, 0) as the body-fixed coordinates which are transformed to the space- 
fixed system by an overall rotation of ip. For a fixed R, these coordinates define a sphere in 
four-dimensions within the following intervals: 

-7r/4 < 9< 7r/4 , 
-7r/2 < (j)< 7i/2 , 

-TT < '0< TT . (39) 

Exchange of two particles is achieved by 

9^-9 , 

and 0^7r/2 — 0, ip ^ ip if 0>O, 

and 0^ -7r/2 - 0, ijj ^n + ip if 0<O. (40) 

With this choice of coordinates, the radial coordinate becomes 

^ r? + r^ (41) 

which is the radius of the sphere in four dimensions. Also, 

X = xiy2 - X2yi = R^ sin(2^)/2 . (42) 

Notice that X depends only on R and 9. Therefore the two-body interaction in the Hamil- 
tonian is independent of the angles and ip. The integrals of motion of the system may be 
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constructed in terms of these new coordinates. The angular momentum operator is given 
by 

d 

which commutes with the Hamiltonian. There exists another constant of motion given by 

d d d d d 

Q = i[x2-^ + y2j; xi- Vit;—] = - ^jTr ■ (44) 

oxi oyi 0x2 oy2 ocp 

Since Q is antisymmetric, acting on a symmetric state produces an antisymmetric state and 
vice versa. We therefore refer to this as a supersymmetry operator (SUSY). The operator 



Q is similar to the SUSY operator discovered in the many-anyon problem by Sen Note 
that the differential operator for both angular momentum and the SUSY operators has a 
very simple form in the hyperspherical coordinates. The states can therefore be labeled 
by the quantum numbers associated with these two operators which we denote by / and 
q respectively. With SUSY, the two-body problem is integrable. (The four constants of 
motion are the Hamiltonian if, the angular part of if, L and Q). Note that we have 
QX = which makes calculations simple. It is easy to check that the bosonic ground state 
of the Hamiltonian has the quantum numbers / and q of the angular momentum and SUSY 
operators equal to zero. 

We would like to emphasize that the eigenstates of the SUSY operator Q are neither 
symmetric nor antisymmetric, unless the eigenvalue q = 0. After finding a simultaneous 
eigenstate of H, L and Q, we can separate it into symmetric (bosonic) and antisymmetric 
(fermionic) parts. These parts are individually eigenstates of but not of Q. Specifically, 
we have Q'^b = q'^F and = q'^B, where B and F denote bosonic and fermionic states 
respectively. Then \E'b ± ^E'f are eigenstates of Q, while and are eigenstates of Q^. 

The two-body Hamiltonian in terms of the hyperspherical coordinates is given by 

1 9^ 3 9 2 
^= -2^d^^RdR-W-^'^^ ^'R^sm\29) ' ^^^^ 

where the operator A^ is the Laplacian on the sphere and is given by 
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cos(2^) de cos2(2e) L dcj^- 



92 2sin(2^) d 1 r 



- + 2sin(2^) 



92 92 



(46) 



The interaction in the Hamiltonian is independent of the angles (j), and depends only on 
R, 9. The operators L and Q commute with the Hamiltonian since they commute with the 
noninteracting [g — 0) Hamiltonian. We thus label the states with the eigenvalues of these 
operators for all gi. Each of these states is four-fold degenerate: Under parity, L —L 
and Q ^ Q and the Hamiltonian is invariant under parity. Therefore the states labeled by 
quantum numbers {I, q) have the same energy as {—l,q). The Hamiltonian is also invariant 
under the transformation fi —fi and f2 ^2- This is a discrete symmetry peculiar to 
this system. Under this transformation L ^ L and Q — > —Q. Therefore the states labeled 
by quantum numbers {I, q) have the same energy as {I, —q). Combining the two we get the 
four- fold degeneracy of the states. Later we will find that the states with (/, q) have the same 
energy as {q, I) since interchanging q and / leaves the differential equation invariant; therefore 
the energy of these two states must be the same. We thus have an eight-fold degeneracy 
for the levels for which |g| and \l\ are nonzero and different from each other. Note that this 
degeneracy is a subset of the degeneracy of the noninteracting system. If |/| = \q\ is nonzero, 
we have a four-fold degeneracy. Finally, there is a four-fold degeneracy between the states 



We are now interested in solving the eigenvalue equation given by H'^ — . Following 
the remarks made in the previous subsection, we may in general write 



The eigenvalue equation separates into angular and radial equations. The angular equation 

is given by 



(±Z,0) and (0,±/) if Z 7^ 0. 



A. Solutions of the eigenvalue equation 



* = F{R) ^0,(f),ip) . 



(47) 




+ 



4^1 



)$ 



/5(/3 + 2) $, 



(48) 



sin2(2^) 
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where (3 > —I, and the radial equation is given by 

The radial equation can be easily solved using the methods outlined in the last section of 



r?| . The solution is given by 

F{R) = R^M{a, b, R^) exp {-R^ /2) , (50) 

where h = 13 + 2 and a = {(3 + 2 — E)/2 and M{a,b, R'^) is the confluent hypergeometric 
function. Demanding that a = —Ur where is an integer, the energy is given by 

E = 2nr + P + 2 . (51) 

Note that (3 is still unknown and has to be obtained by solving the angular equation. Never- 
theless the tower structure of the eigenvalues built on radial excitation of the ground states 
is obvious from the above. 

The angular equation (^8[) may be solved with the ansatz 

= P{x) exp(ig0) exp{il^) , (52) 

where x = sin(2^) and /, q are the state labels in terms of the integer valued eigenvalues 
of the angular momentum and SUSY operators. The angular equation then reduces to a 
differential equation in a single variable x for the function P{x): 

^,£P dP 1 , , , Qi 0(0 + 2) 

1 -x^)—-- 2x— - — + 2xql + l^]P-^P+ t^^ti^p = q . 53 

dx^ dx 4:{1~ x^) ^ x^ 4 

Note that the equation has four regular singularities at x = 0, 1, —1, 00 (the singularity at 
00 does not play any role since x is bounded). Therefore the solution is of the form 

P{x) = |xP(l-x)^(l + x)"e"'^'^(a;) . (54) 

One can now fix a, b, c to cancel the singularities. We find that 

b = and c = — . (55) 

4 4 ^ ^ 
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Since / and q are integer valued, the values of h and c are restricted. The other exponent a 
is given by 

a(a — 1) = gi = g{g — 1) with a > , (56) 

where we have already defined gi through Eq. ([^ in terms of g. Note that we have used the 
symbol a instead of g. As shown in the Appendix, we have to take a = g if g > 1/2. But if 
g < 1/2, we have to generally consider a linear superposition of solutions with a equal to g 
and 1 — g (more on this later). 

We finally arrive at the required differential equation from which the eigenvalues are 
determined, 

(1 - x^)— r + 2[a/x - (6 - c) - (a + 6 + c + l)x] — 
dx'^ ax 

+ [ ' - (a + 6 + c + 1/2)2 ^ 2a(c - h)/x ] O = . (57) 

For (? = 0, the solutions are simply Jacobi polynomials and the full solution for the angular 
part is given in terms of the spherical harmonics on a four-dimensional sphere. In general, 
this differential equation is known as the Heun equation whose solutions Q"-'^''^{x) are char- 
acterized by the so-called P-symbols The Heun equation is exactly solvable if either I 



or q vanishes, i.e., if 6 = c as discussed in the next subsection. The equation is also exactly 
solvable at an infinite number of isolated points in the space of parameters (a,6, c). These 
are isolated points because if we vary a slightly away from any one of them, the equation is 
not exactly solvable. Note that 6, c take discrete values and cannot be varied continuously. 

B. Polynomial solutions 

Let us first consider a class of solutions which are polynomials in x. We may then write 

e(x) = Ckx'^ ' (58) 

fc=0 

where we may define Cq = 1. Substituting this in the differential equation for 6, we see 
that the C{.s satisfy a three-term recursion relation given by 



{k + 2){k + 1 + 2a) Ck+2 - 2(6 - c){k + 1 + a) Ct+i 

+ [^-^^^^-{a + h + c + k + l/2f]Ck = 0, (59) 

which is in general difficuh to solve. However there are two special cases when polynomial 
solutions are possible, (i) For b = c, this reduces to a two-term recursion relation which 
can be easily solved to obtain all the energy levels. This is an example of a Conditionally 
Exactly Solvable (CES) problem ||21]] in which the full spectrum is exactly solvable for some 
special condition (like b = c here), (ii) The other case is when the coefficient of Ck is zero 
with k = p (where p > 1), i.e. 

E = 2nr + 2a + 2b + 2c + 2p + 2 , (60) 

in which case one has a polynomial solution of degree p. This is an example of a Quasi- 
Exactly Solvable (QES) problem when only a finite number of states are exactly solvable 
for some given values of the parameters. As far as we are aware, this is the first example 
where both CES and QES solutions exist in the same problem. We now discuss both types 
of solutions in detail. 

(i) CES-type Solutions: To see the solutions explicitly, define 

y = x"^ . 

In terms of the variable y the differential equation is written as, 

yil-y)^ + 2[a'-ib' + c' + l)y]^-b'c'Q = 0, (61) 

where 

a = a + 1/2 , 

b' = i[a + 26 + l + |] , 

c' = i[a + 26-|]. (62) 

This is now a hypergeometric equation whose solutions are given by F{b',c',a';y). Note 
that we need to concentrate only on the solutions for < y < 1. The hypergeometric 
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series terminates whenever 6' or c' is a negative integer or zero. In our case h' is always 
positive hence the solutions are given when c' = — m, where m is an integer. Therefore 
j3 = 2m + 2a + 46. The energy eigenvalues are obtained by substituting this in Eq. (pTf), 



E = 2nr + 2a + Ab + 2m + 2 . (63) 

Since 6 = c, we must have either / = or g = 0. The energy varies linearly with a as in 
the case of the exact solutions of the many-body problem. The expression for a in terms 
of g will be clarified in the next subsection. It might be tempting to conclude that all 
polynomial solutions vary linearly with g. In fact that is not so as can be seen from the 
following examples. 

(ii) QES-type Solutions: We can have polynomial solutions of degree p (where p > 1) if 6, c 
and g satisfy some specific relations. Let us consider the case 

Q{x) = l + 6x . (64) 

This is a solution if 

"^^ 1 (65) 



{b - C)2 - 1 

and 6 = b — c. We have implicitly assumed that b is not equal to c since otherwise this 
solution is trivial and of CES type. (Once again, a can be equal to either g or 1 — g as 
discussed in the previous subsection). Then /3 = 2a + 26 + 2c + 2, and the full spectrum is 
given by 

E = 2nr + 2a + 2b + 2c + 4 . (66) 

We should point out that a solution of this kind is only possible for fairly large values of / 
and q; the minimum values needed are |/| = |g| = 3, in which case a = 1/5. 
Similarly, there is a polynomial solution of the form 

e{x) = l + Sx + ex^ , (67) 

if 
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3 / , 1 
a = 2a - - + \ Aa-' - a + - , 
2 V 4 

where a = ^ — - . (68) 

The spectrum is given by 

E = 2nr + 2a + 2b + 2c + 6 . (69) 

These expressions for the energies are nonhnear in a because of the constraints (|65D or 
(|68|). We should however caution that these are isolated solutions since b and c can only 
take discrete values; hence the above solutions do not vary smoothly with a. In general, 
solutions similar to the above may be constructed for every degree p of the polynomial; 
the corresponding energies are given by Eq. (|60|) where a is given by a function of b and c 
which can be derived by solving p+1 recursion relations obtained from Eqs. ( |59D by setting 
fc = -l,0,l,...,p-l. 



C. Numerical analysis 

We now consider the numerical solution of some low-lying states of the two-body problem 
since the polynomial solutions described in the previous subsection do not exhaust the full 
spectrum. 

The noninteracting limit of the system is g = where we have the solutions corresponding 
to a four- dimensional oscillator. These are simply the spherical harmonics on a four-sphere 
Yk^i^q, where k = 0, 1,2, ■ ■ ■ and |/|, |g| < k label the states. When the degeneracy of these 
states is taken into account, all the states in the noninteracting limit g = are completely 
specified. We now demand that the wave functions and energy levels should vary continu- 
ously with the parameter g which is the interaction strength. We also require that the wave 
functions should not diverge at any value of x in the interval [—1,1]. 

For given values of 6 7^ c, we can numerically find the energy levels in two different ways. 
We can diagonalize the differential operator in ( ^31) in the basis of the noninteracting {g = 0) 
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states, or we can solve the differential equation (|5^) or (^7^ directly for each state. We have 
used both methods and will present the results below. 

In order to proceed further, it is necessary to clarify the dependence of a on g in Eq. 
(^). By the arguments given in the Appendix, a = g if g > 1/2. If g < 1/2, a can be equal 
to either g or 1— g; for any given state, we choose the value which is continuously connected 
to the noninteracting solution at g = 0. Namely, any solution of the above kind must, at 
g = 0, go either as 1 or as x near x = 0; we then choose a = (yfora = l — (7in the two cases 
respectively for < < 1/2. In particular, an exact solution which, near x = 0, goes as 1 
at g = will go as \x\^ for all g > 0. An exact solution which goes as x at g = will go 
as upto g = 1/2 and then as \x\^ for g > 1/2. As a result solutions of this form are 

discontinuous at g = 1/2. In general, when solving numerically for the non-exact solutions, 
we have to allow a superposition of both \x\^ and For g > 1/2, however, a must be 

equal to g. 

When solving the differential equations ( pSf ) or (|57D, we have to consider the two regions 
separately. According to the rules discussed in the Appendix, for < g < 1/2, we take the 
function P{x) to go as 

P(x) = + d sgn(x) , (70) 

near a; = 0, and we vary both the coefficient d as well as j3 in Eq. (|53|) till we find that the 
solution of (|5^) does not diverge at x = ±1. Both d and P depend on g; it is possible to 
determine the limiting values d{0) and d{l/2) as follows. At g = 0, suppose that the Jacobi 
polynomial is normalized so that P = 1 and P' = Cj at x = 0. On the other hand, from Eq. 
(^), we see that Ci = b — c for any nonzero g, no matter how small. By taking the limit 
g ^Oin Eq. (^, we therefore find that d{0) =Cj-b + c. At the other end, d{l/2) = ±1 
because, as we will show next, P{x) must vanish for either x > or for x < for all (7 > 1/2 
(except for the exact polynomial solutions discussed previously). 

For g > 1/2, we must take the solution of the differential equation to go either as (a) 
P(x) = for X < and ~ x^ for x small and positive, or as (b) P(x) = for x > and 
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~ (— x)^ for X small and negative. In either case, we vary the energy till we find that the 
solution of ( |57D (with a = g) does not diverge at x = 1 and —1 for (a) and (b) respectively. 
It is interesting to note that for each such solution, P{x) vanishes identically in one of the 
half-intervals [—1,0] or [0, 1]. Note also that if B"'^''^(x) is a solution which is only nonzero 
for X > 0, then G'^''^'^(x) = Q"''^''^{—x) will be a solution which is only nonzero for x < 0; 
further, the two solutions will have the same energy. 

In contrast to the above method for solving the Heun equation ( ^7|) directly, the nu- 
merical diagonalization procedure for finding the eigenvalues involves solving the eigenvalue 
equation (^81) . The basis for diagonalization is provided by the eigenstates of the Laplacian 

on S^, namely the spherical harmonics on a four-sphere Yf^^i^g. For non-zero interaction 
strengths, the singular interaction is handled by multiplying the non-interacting eigenstates 
by |xp. The resulting basis is nonorthogonal, and the diagonalization procedure is fairly 
straightforward though cumbersome. We have truncated the basis such that the highest 
energy state in the basis is lOO^a;. The results for the low lying states in the spectrum 
obtained through both methods are displayed in Fig. 1 (for < g < l)and Fig. 2 (for 
g > 1). For g < 1/2, it is more convenient to use the diagonalization procedure since the 
direct solution of Eq. ( pB] ) requires one to numerically fix two separate parameters d and (3. 
On the other hand, it is easier to solve the Heun equation ( ^7|) for g > 1/2 since one has to 
fix only one parameter (3. In general we have used both methods to arrive at the spectrum 
of low-lying states. For small values of g (< 1), however, the solutions of the differential 
equation produce eigenvalues which are somewhat smaller than the ones obtained by the 
diagonalization procedure. For large values of g, however, there is no perceptible difference 
between the results from the two methods. 

Fig. 1 shows the energies for some values of I and q. Each level is labeled by {l,q]D), 
where D is the degeneracy of the level away from the noninteracting limit; the degeneracy is 
computed by counting the allowed values of ±/ and ±g using the parity and supersymmetry 
transformations for a given level. A subscript 'r' on the label (/, q, D) denotes the radial 
excitation which is simply inferred from the existence of the towers. The bosonic ground 
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state has the predicted behavior for all (7; it is linear with a slope 2 as a function of g. The 
corresponding wave function goes as \x\^ as x — 0. In contrast, the level (0,0; 1) starting 
at -E = 4 has an entirely different behavior . It is exactly solvable for all g. According to 
the previous discussion, dE/dg = —2 for (7 < 1/2 (with the wave function going as 
for small x) and dE/dg = 2 for g > 1/2 (with the wave function going as thus dE/dg 
is discontinuous at g = 1/2. 

At g = 0, the slopes dE/dg for all the levels can be calculated using first-order per- 
turbation theory as shown in the next subsection. For large values g, we find that all the 
energies converge to 2g plus even integers as shown in Fig. 2. This amazing behavior can 
be understood using the WKB method as shown in subsection E. 



D. Perturbation theory around g = 

It is interesting to use perturbation theory to calculate the changes in energy from g = 0, 
and to compare the results with the numerical analysis. We will only describe first order 
perturbation theory here, and the example we will consider is the fermionic ground state 
which is doubly degenerate for N = 2. 

In general, naive perturbation theory fails at g = because most g = eigenstates do 
not vanish as Xij 0; hence the expectation value of l/Xfj diverges. This problem can be 
tackled by using a special kind of perturbation theory first devised for anyons [^. We will 
first describe the idea for particles and then specialize to = 2. Instead of solving the 
equation = E'^, where H is given in Eq. (^, we perform a similarity transformation to 
H = X^^HXlj and § = X^^^, where 

Xn = W \Xij\^ . 

i<j 

We then find that H = Hq + V, where Hq is the noninteracting Hamiltonian (with gi = 
g2 = 0), and 
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The first order changes in energy may now be obtained by calculating the expectation values 
(or matrix elements, in the case of degenerate states) of V in the zeroth order (noninteract- 
ing) eigenstates. These expectation values can be shown to be convergent for all states. Note 
that V only contains two-body terms. Although V is not hermitian, it is guaranteed that 
its expectation values are real because the original problem has a hermitian Hamiltonian H. 
For = 2, we find that 

o . 1 ^ cot 20 d , 

in hyperspherical coordinates. Note that V commutes with both L and Q, so that we only 
have to consider its matrix elements within a particular block labeled by the eigenvalues I 
and q. Let us now use ( |72D to compute the first order change in the states which have E = 3 
at g = 0. There are four such states, with / = ±1 and g = ±1 (labeled (1, 1;4) in Fig. 1); 
two of these states are actually the ground states of the two-fermion system. Due to parity 
and SUSY, these four states remain degenerate for all g. Hence it is sufficient to calculate 
the first order change in the state with, say, {l,q) = (1,1)- Since this state is unique at 
g = 0, we only need to do non-degenerate perturbation theory with V. The normalized 
wave function for this state is 

^ = (cos - sin 0) exp[i(0 + ^)] /?exp[-i?V2] . (73) 

vrV 2 

We now obtain the expectation value 

/ dR cos(20) d9 d(p di) m* = g . (74) 

Jo J-Tt/A J--K/2 J-TT 



We can see from Fig. 1 that this gives the correct first order expression for the energy 
E = 3 + g near g = for the states labeled (1, 1;4); their first radial excitations (1, l;4)r 
therefore have E = 5 + g. 

We can similarly calculate the first order expressions for the energies near g = for all 
the other levels shown in Fig. 1. We find that E = 2 + 2g for the lower state labeled (0, 0; 1) 
(i.e. the bosonic ground state); E = 4 + 2g for its radial excitation (0, 0; 1)^., and the states 
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(2, 0; 4); £^ = 4 - 2g for the upper state (0, 0; 1); E = 4 for the states (2, 2; 4); E = 5 - ^/2 
for the states (1, 1; 4) and (3, 3; 4); and E = b + 3(7/2 for the states (3, 1; 8). 



E. Large-(7 perturbation theory 

We can study the solutions of Eq. (|57D for large values of g by using an expansion in 
1/ g. For any value of h and c, we will only study the lowest energy and we will calculate 
the leading order terms in E and the wave function G(x). We first note that the terms of 
order g"^ in ([57| ) can be satisfied only ii E = 2g + 0{1). Next, we assume that E and have 



WKB expansions of the form 



E = 2^ + 26 + 2c + 2 + /o + - + 0{llg^) , 

9 

e = exp[^i;o(a;) + ] • (75) 



The boundary condition 6 = 1 implies that Wo{0) = Wi{0) = 0. To order g, Eq. (|57|) gives 
the first order differential equation 



2n dwo ^ fo 



(l-x^)^ = b-c-'-^x. (76) 

We now look at solutions which are nonzero only for x > 0. We demand that G should 
neither diverge nor vanish (since the lowest energy solution should be node less) anywhere 
in the range < x < 1. Hence the functions Wq and Wi should not diverge to oo or —oo in 
that range. This fixes fo = 2(6 — c), so that 

E = 2^ + 46 + 2 = 2^ + + + 2 , 
and e = (1 + xf-^ . (77) 

Similarly, there are solutions which are nonzero only for x < 0. These have 

E = 2^ + 4c + 2 = 2^ + |/ - g| + 2 , 
and e = (1 - xy-^ . (78) 
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We now go to order 1 in Eq. (|57D. For solutions which are nonzero only for x > 0, we find 
that 

E = 2g + Ah + 2 + '—— = 2^ + |/ + g| + 2 - , 

9 4:9 

and e = (1 + xf-' exp[ + ^ + (ln(l + x) - ] . (79) 

-L I X 

We can similarly find solutions which are nonzero only for x < 0, by changing x — > —x and 
interchanging b ^ c, i.e. I + q I — q and Iq —Iq, in Eq. ([79|) . 

We see from Fig. 2 that these formulae correctly describe the leading behavior of E. 
In fact, the laige-g behavior is already visible in Fig. 1, for some states, as we approach 
g = 1. The various levels shown in that figure have the following WKB energies; E = 2g + 2 
for both the (0, 0; 1) states (one of these is the bosonic ground state and the other is the 
fermionic ground state for > 1/2 as discussed later); E = 2(7 + 2 + 1/4(7 for the states 
(1, 1; 4); E = 2^ + 2 + 1/^ for the states (2, 2; 4); E = 2^ + 2 + 9/Ag for the states (3, 3; 4); 
E = 2(7 + 4 — 1/4(7 for the states (1, 1; 4); E = 2g + A for the radial excitation (0, 0; 1)^ and the 
states (2,0; 4); E = 2g + A + l/Ag for the radial excitations (1, l;4)r; and E = 2g + A + 3/4:g 
for the states (3, 1; 8). We have also checked that the leading order wave functions in Eqs. 
([f9|) agree remarkably well with the correct wave functions 0(x) obtained by solving the 
Heun equation ([57|) even if g is not very large. 

It is easy to see from Eqs. ([77| - [78D that for large g, the ground state and also the excited 



states become infinitely degenerate. This is so because one can choose the quantum numbers 
/ and q in infinitely many ways such that the energies are the same as (7 — oo. Further, the 
spacings now become twice the spacing at (7 = since / and q have the same parity mod 2. 

The large-(7 behavior therefore displays a remarkable similarity to the problem of a par- 
ticle in an uniform magnetic field where the Landau level spacing is twice the cyclotron 
frequency, and each level is infinitely degenerate. 
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F. Fermionic Ground State Energy 



The fermionic ground state energy has a very unusual behavior as can be seen from Fig. 
1. For < g < 0.367, the ground state energy monotonically and nonhnearly increases from 
3 to 3.266 along the curve labeled (1,1; 4). Beyond this point, for 0.367 < g < 0.5, the 
ground state energy monotonically and linearly decreases from 3.266 to 3 along the upper 
curve (0,0; 1) satisfying E = 4 — 2g. For g > 1/2, the fermionic and bosonic ground state 
energies are identical and are given by the curve (0,0; 1) which satisfies Eq = 2 + 2g, i.e., 
both the ground states monotonically increase with g. Thus the fermionic ground state 
consists of three pieces as a function of g, while the bosonic ground state is given by the 
single line Eq = 2 + 2g for all g > 0. 

For two particles, one can understand why the fermionic and bosonic ground state en- 
ergies are identical foTg> 1/2 as follows. In this range of g, the ultraviolet potential near 
x = is infinite and it prevents tunneling between the regions x > and x < (see the 
Appendix). For two identical particles, an exchange necessarily takes us from a region with 
X > to a region with x < according to (^Ol). If tunneling between the two regions is 
forbidden, it becomes impossible to compare the phases of the wave function of a given con- 
figuration of the two particles and the wave function of the exchanged configuration. Thus 
it is impossible to distinguish bosons from fermions if g > 1/2, and their energy levels must 
be identical. 

It is possible that the same argument will go through for more than two particles; however 
we need to understand the ultraviolet regularization of the three-body interactions properly 
in order to prove that rigorously. If the argument holds, then we would have the interesting 
result that the iV-fermion ground state energy is also given by (H) for g > 1/2, while it may 
show one or more level crossings for g < 1/2. 
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IV. DISCUSSION AND SUMMARY 



To summarize, we have studied a two-dimensional Hamiltonian whose eigenstates have 
a novel two-particle correlation. We have shown the existence of several classes of exact 
solutions in the many-body problem. We have analyzed the two-particle problem in detail 
and shown that it is completely solvable by reducing it to an ordinary differential equation 
in one variable which can be solved exactly for a subset of states and numerically otherwise. 
The two-body problem is integrable since there are four constants of motion in involution. 
We have also discussed perturbation theory for both small and large coupling strengths. In 
the strong interaction limit, the system simplifies and bears a remarkable resemblance to 
the Landau level structure. 

We have also clarified in the Appendix the ultraviolet prescription which is required to 
make sense of an inverse- square (singular) potential especially at small coupling strengths. 
In particular, we emphasize that it is in general not sufficient to specify that the wave 
functions are regular and square integrable to obtain an energy spectrum uniquely when 
dealing with singular interactions. In some domains of the coupling strength, we also need 
to specify the ultraviolet regularization to make complete sense of the results. We do this 
by demanding that as the parameter (7 — > 0, the energy levels should smoothly approach the 
known noninteracting levels. We believe that this discussion is quite general and may have 
a wider applicability to Hamiltonians with singular interactions. 

Interesting problems for the future would be to extend this analysis to more than two 
particles, and to find an application of our model to some physical system. Recently, we 
have come to know that our model has been generalized to three (and higher) dimensions 



with novel three-body (and many-body) correlations p4 
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Appendix 

We begin directly from Eq. (pS]). Given the real number g > satisfying gi = g{g — 1), 
P{x) could go, as |a;| —>■ 0, as either \x\^ or |a;|^~^ or even as a general superposition of the 
two powers. We therefore need to define the problem more carefully in order to pick out a 
desired solution . 



As mentioned above in the text, we demand the following. Firstly, the limit g = should 
give all the noninteracting two-particle solutions, both bosonic and fermionic. Secondly, all 
the wave functions and energies E should be continuous functions of g, but the first derivative 
of E need not be continuous (indeed dE/dg is not always continuous at (7 = 1/2 as we saw 
earlier). Finally, for (7 > 1, the wave function should go as and not as which 
diverges at x = 0. 

iFiom. these three requirements, it is clear that for g > 1/2, the wave functions must go 
purely as \x\^, whereas for g < 1/2, the wave function could go either as \x\^ or or a 

superposition of the two. 

We will now show that we can satisfy the above requirements if we redefine the problem 
with a different potential in an ultraviolet region |x| < Xq- We take the potential to be 

V(x) = ^ for |a;| > Xo , 

X'^ 



xl 



for \x\ < X, 



o 5 



where u tanhw = g if < g < 1/2 , 

and M = cx) if ^ > 1/2 . (80) 

Eventually, of course, we have to take the limit Xq — ^ to recover our original problem. 
Note that the potential in the ultraviolet region is not symmetric under g ^ 1— g for g < 1. 
Hence the energy spectrum does not have this symmetry. 

To see why Eqs. ( pUD work, we note that the wave function, for |x| slightly greater than 
Xo (where Xo is much smaller than any physical length scales like the width of the harmonic 
oscillator potential), is generally given by 
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P(x) = + d+ ^ if x> Xo 
and P{x) = {-xY + i-xY'^ if x< -Xo . (81) 

(For the exceptional case = 1/2, we have to replace \x\^ and by \x\^^'^ and \x\^^'^ In 

respectively). 

Now consider the first case in (POD, i.e., < (7 < 1/2. Since the energy E is much less 
than the potential in the inside region |x| < Xo (this is necessarily true for any finite value 
of -E as So ^ 0), the wave function in that region is given by 

11 

P{x) ~ cosh [ ( — + 0{xo) ) X + 6] , (82) 

Xq 

where 6 can be a complex number, and the term of 0{xo) arises from the energy E which is 
much less than {u/xoY- We now match the wave function and its first derivative, or, more 
simply, the ratio P'{x) / P{x) aX x = Xq^ e and aX x = —Xq ± e, where e is an infinitesimal 
number. We then find three possibilities. 

(i) The wave function may be even about x = 0. Then 5 = 0, and dj^ = must vanish as 

0. (The behavior of d± can be deduced by equating the terms of 0{Xq ^) and 
0{xo) in P'/P at X = Xo±). In the limit Xo — > 0, therefore, the wave function goes purely 
as \x\^. 

(ii) The wave function may be odd about x = 0. Then 6 = i7c/2, and c?_|_ = d^ must diverge 

— >■ 0. The wave function is proportional to sgn(x) ^ in that limit. 

(iii) In the general asymmetric case, we find that we must have 6 of order x^"^^, and d+ = 
— (i_ = d of 0(1). (This is found by equating terms of 0(a;^^) and 0{x~'^^) in P'/P at 
X = ±Xo). The wave function is therefore a superposition of the form 

P(x) = \x\^ + rfsgn(x) . (83) 

The cases (i) and (ii) arise if either / or g is zero in Eq. (0), since the equation is 
invariant under x —x in that case. This is precisely when b = c and the equation is 
exactly solvable. We thus see that the even solutions go as \x\^, while the odd solutions go 
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as If neither / nor q is zero, i.e. b ^ c, we have case (iii) where a superposition of the 

two powers are required. 

The second case in (pO|), i.e. g > 1/2, is relatively simpler to analyze since the wave 
function must be zero in the inside region |a;| < Xo- On imposing this condition on the wave 
function in the outside region, we see that both and (i_ must vanish as Xq — > 0. Hence 
the wave function will go purely as \x\^ in that limit. However, since there is no tunneling 
possible through the infinite barrier separating x > Xo from x < —Xo, we will generally have 
wave functions which are nonzero only for x > Xo or only for x < —Xo- This is indeed true 
as we saw earlier for the solutions of the Heun equation for g > 1/2. 

We would like to emphasize that the relation wtanhw = g in Eq. (|80D is absolutely 
essential in order to have the possibility of P{x) ~ \x\^ for g < 1/2. li u were to take any 
other value, we would find that P{x) necessarily goes as in the limit Xq — 0. A similar 

fine tuning of u is necessary in the GSM for < 1/2. Incidentally, the strongly repulsive 
potential in the ultraviolet region explains the peculiar result that the bosonic ground state 
energy increases monotonically with g even though the potential away from the ultraviolet 
region becomes more and more attractive as g goes from to 1/2. One can show from Eq. 
(^) that the integrated potential j\ ldxV{x) is actually positive and large if Xo is small, 
and it increases as g varies from to 1/2. 

Several comments are in order at this stage, 
(i) A similar fine tuning is also required in the CSM ii g < 1/2 is to be allowed as has been 



done by several people |T3[. Historically, both Calogero |]T| and Sutherland restricted 
themselves to (7 > 1/2. In a sense they could do that since the free fermion limit corresponds 
to g = 1; thereby they avoided the problem with g < 1/2. However, one cannot reach the free 
bosonic limit smoothly in that case. Both these authors believed (7 < 1/2 to be unphysical 
because they chose a particular regularization. What we have argued here is that one can 
choose an alternative regularization (called the resonance condition by Sutherland) which 
allows one to go continuously all the way upto g = and hence reach the free bosonic 
limit continuously. We have not seen this being clearly stated in the CSM literature before, 
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although Scarf |]T2| discusses this issue in a different problem containing the inverse-square 
potential. 

It is worth noting that CSM has only two-body interactions. Henca the entire discussion 
here is also valid in the many-body case in CSM. This is in contrast to our problem where, for 
N > 2, one also has to analyze the ultraviolet regularization of the three-body interactions. 

(ii) One important consequence of our regularization is that many of the states have dis- 
continuities in dE/dg aX g = 1/2; the fermionic ground state also has a level crossing at 
g = 0.367. Further, for each value of gi < 0, there are two possible ground states since the 
ultraviolet regularization depends on g and not on gi. 

(iii) Put differently, our Hamiltonian H has several self-adjoint extensions (SAE) for each 
value of gi. What we have done is to choose a particular SAE for gi > and two different 
SAE for gi < 0. As a result, we have found that for every value of gi in the range —1/4 < 
gi < 0, there are two possible ground state energies since, as seen above, the SAE depends 
on g rather than on gi. Actually, there is an even more general SAE possible for any value 
of gi where another real parameter (besides g) has to be introduced; however we shall not 
discuss that here. 
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FIGURES 

FIG. 1. Spectrum of some low energy states in the two-body problem as a function of the 

interaction parameter g. The states arc labeled by {l,q;D), where the first two entries denote the 
two angular quantum numbers and the last entry shows the degeneracy of the level. Further the 
radial excitations are denoted by the subscript 'r'. 

FIG. 2. Spectrum of states as a function of g showing the strong coupling behavior. We have 
shown E — 2g instead of the energy itself. 



36 




Fig.l of Bhaduri, Khare, Law, Murthy & Sen. 




Fig.2 of Bhaduri, Khare, Law, Murthy & Sen. 



